An artificial neural network application 
on nuclear charge radii 



Serkan Akkoyun, Tuncay Bayram, S. Okan Kara and Alper Sinan 

Department of Physics, Cumhuriyet University, Sivas, Turkey 
Department of Physics, Sinop University, 57000 Sinop, Turkey 

Department of Physics, Nigde University, Nigde, Turkey 
Department of Statistics, Sinop University, 57000 Sinop, Turkey 



Abstract. The artificial neural networks (ANNs) have emerged with 
successful applications in nuclear physics as well as in many fields of 
science in recent years. In this paper, by using (ANNs), we have con- 
structed a formula for the nuclear charge radii. Statistical modeling of 
nuclear charge radii by using ANNs has been seen as to be successful. 
Also, the charge radii, binding energies and two-neutron separation en- 
ergies of Sn isotopes have been calculated by implementing of the new 
formula in Ifartree-Fock-Bogoliubov (HFB) calculations. The results of 
the study shows that the new formula is useful for describing nuclear 
charge radii. 
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1 Introduction 

Nuclear charge radii and, more generally, distributions of charge-density give 
direct information on the Coulomb energy of nuclei. Because of this reason, it 
has gained attention for the nuclear mass formulas in recent six decades |T|. As 
is well known, the nuclear charge radii is a fundamental property of atomic nu- 
clei [213] . It can be measured by various methods based on the electromagnetic 
interaction between the nucleus and electrons or muons. Widely used methods 
are measurements of transition energies in muonic atoms, elastic electron scat- 
tering experiments, -fCaX-ray and optical isotope shifts. A detailed discussions 
on these techniques for measuring of the rms (root-mean-square) charge radii of 
nuclei can be found in [1]. The latest advances in experimental techniques have 
provided accessibility to experimental nuclear charge of more nuclei far from 
/3-stability line. 

The nuclear charge radii related with exotic phenomena such as skin and halo 
has become a hot topic [5]. With its accuracy, the study of the nuclear charge 
radii is very important for a better understanding of not only the proton distri- 
bution in nuclei but also the skin and halo. For these reasons, if one get a simple 
and reliable formula for nuclear charge radii, it can be provides information for 
exotic nuclei and the effective nucleon-nucleon interaction. For these purpose, a 
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number of nuclear charge formulae have been proposed as to be based on A-^^^ 
and Z^^^ dependence. An excellent comprehensive study on nuclear charge for- 
mula obtained from fitting of experimental data can be found in 5j and reference 
therein. However, in these studies experimental data were limited. New exper- 
imental data over more than 800 nuclei have been updated, recently Using 
the new data, the parameters of the nuclear charge formulae can be revised. In 
the present work, we used feed-forward artificial neural networks (ANNs) as a 
different approach for obtaining a reliable and simple nuclear charge formula. 
As is well known, feed-forward ANNs may be viewed as an universal non-linear 
function approximator In recent years, ANNs have been used in many fields 
in nuclear physics as in the other fields, such as determination of one and two 
proton separation energies [7], developing nuclear mass systematics [8|, identi- 
fication of impact parameter in heavy-ion collisions [9)10|llj . estimating beta 
decay half-lives [H] and obtaining potential energy curves [T3]. The fundamen- 
tal task of the ANNs is to give outputs in consequence of the computation of the 
inputs. We have obtained new mass dependent simple nuclear charge formula by 
using all the data from ANN whose inputs were taken from the latest experimen- 
tal data. The obtained new formula is presented and used in HFB calculations 
for Sn nuclei. 

The paper is organised as follows. In Section 2, ANNs formalism and Hartree- 
Fock-Bogoliubov (HFB) method are given briefly. In Section 3, the results of the 
study and discussions are presented. Finally, a summary is given in Section 4. 

2 Theoretical Framework 
2.1 ANNs 

ANNs are mathematical models that mimic the human brain. They consist of 
several neurons which are processing units. The neurons are connected each 
other via adaptive synaptic weights [2]. The ANN is composed of three main 
layers. The first layer corresponds to input layer, the intermediate layer is called 
hidden layer and the last one is output layer. In this study, one input layer with 
a neuron, one hidden layer with forty four neurons (h=44) and one output layer 
with a neuron ANN was used. The total adjustable weights was 2h=88. 

The input layer neuron receives the data from outside and the output layer 
neuron gives the results. The data is transmitted via weighted connections be- 
tween the neurons. The tangent hyperbolic (sigmoidlike) and linear functions 
were used for hidden and output units respectively. It has been proven that one 
hidden layer and sigmoidlike activation function in this layer are sufficient to 
approximate any continuous real function |15j . 

The use of ANN is as two-step process, training and test stages. In this work, 
a back-propagation algorithm with Levenberg-Marquardt jl6ll7j for the training 
of the ANN was used. In the training stage, ANN modifies its weights until an 
acceptable error level between desired and predicted outputs is attained. The 
error function which measures this difference was mean square error (MSE). 
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After an acceptable error level, the trained ANN is tested over the data of 
interest. For further and general background for the ANN, the reader is referred 

to [Tg. 

2.2 HFB Method 

In the HFB formalism, a two-body Hamiltonian of a system of fermions can 
be interpreted in terms of a set of annihilation and creation operators. The 
ground-state wave function is described as the quasi-particle vacuum and the 
linear Bogoliubov transformation provides connection between the quasiparticle 
operators and the original particle operators. The basic building blocks of the 
HFB method are the density matrix and the pairing tensor, and expectation 
value of the HFB Hamiltonian could be expressed as an energy functional (De- 
tails can be found in |2ll8j l. In terms of Skyrme forces, the HFB energy has 
the form of local energy density functional contains the sum of the mean field 
and pairing energy densities. These fields can be calculated in the coordinate 
space |2|18j . In this work, HFB equations have been solved by expanding quasi- 
particle wave functions that conserve axial symmetry and parity on a harmonic 
oscillator basis expressed in coordinate space proposed by Stoitsov et al. |18j . 
For pairing, Lipkin-Nogami method is implemented by performing the HFB cal- 
culations with an additional term included in the HF Hamiltonian. To obtain 
ground-state properties of Sn isotopes prescription of Ref. 18J was taken into 
account. A number of effective Skyrme forces can be found in literature for cor- 
rect prediction of the nuclear ground-state properties. In this study, the Skyrme 
force SLy4 [19^ has been used. 

3 Results and Discussions 

In this study, nuclear masses and experimental nuclear radii of ^900 nuclei have 
been used in ANNs for fitting of the charge radii. The input of the ANN was 
atomic masses of the nuclei and the output was nuclear radii. The whole data 
was partitioned into two seperated sets, %80 for the training stage and the rest 
for the test stage. After the training of the network, it was tested over the test 
set data which have never seen before by the network. As can be seen in the 
figure [T] the ANN predictions agree exceptionally well with experimental values. 

This obviously indicates that the test set ANNs have consistently generalized 
the training set fittings. The MSE values were about 0.0027 and 0.0029 for the 
training and test stages respectively. After testing of the trained network, the 
whole data from ANN outputs was used for employing least square fitting. It 
should be noted that mass dependent form of nuclear charge radii given by the 
formula Rc — r^A^ has been considered for fitting. The novelly obtained formula 
is given by 

Rc = 1.234A°-28 (1) 

where Rc and A is the rms charge radii and mass number of nuclei, respectively. 
This formula are in good agrement with the experimental data. The root mean 
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Fig. 1. The experimental (Exp.) and ANN test set predictions for nuclear radii for 
several isotopes ranging from A=6 to 240. 
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square deviation between charge radii obtained from the formula and experi- 
mental ones is a 1.42:10^^. Beside, based on the formula we calculated the 
ground-state properties of nuclei such as binding energy per nucleon, two- neutron 
separation energy and nuclear charge radii in HFB model for Sn isotopes. 
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Fig. 2. Binding energies per nucleon (upper panel) and two-neutron separation energies 
(lower panel) for Sn Isotopes calculated with the spherical HFB code are compared with 
experimental data (Exp.) (from Ref. [20]). 



The calculated binding energies for even-even '^^""-'^ Sn isotopes within the 
HFB method with the Skyrme force SLy4 are shown in the upper panel of 
Figure |2] Also, available experimental binding energies per nucleon taken from 
Ref. [20] are shown for comparison. As can be seen in the upper panel of the 
Figure [2j the predictions of HFB method are in good agrement with the ex- 
perimental data. Maximum difference between the calculated and experimental 
values is only ~0.04 MeV at neutron number N = 50. In the lower panel of 
the Figure [2j the calculated two- neutron separation energies of Sn isotopes and 
experimental ones are shown. Abrupt decrease of the two-neutron separation 
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energy at neutron number N — 82 indicates that this nuclei has a shell closure. 
Thus, shell effect is clearly visible in our calculations. 
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Fig. 3. The calculated rms nuclear charge radii within the HFB method with the 
Skyrme force SLy4 for Sn isotopes obtained by implementing the new formula in HFB 
code. Also, the experimental rms charge radii of Sn isotopes [4] are shown for compar- 
ison. 



The calculated rms nuclear charge radii of Sn in HFB method by 

implementing the new formula and available experimental charge radii of Sn 
isotopes taken from Ref. [3] are shown in Figure|3] As can be seen in the Figurejsj 
the HFB method with the SLy4 parameters reproduced rms charge radii of Sn 
isotopes well. It should be noted, however, that ongoing from neutron number 
A'^ = 58 to neutron drip line, the difference between calculated and experimental 
charge radii for Sn isotopes is rising. However, maximum difference is only 0.032 
fm at neutron number N — 82. 



4 Summary 

ANNs have been applied on new experimental nuclear charge radii data to obtain 
simple charge radii formula based on mass dependence. The formula is founded 
as to be successful. By using the formula in HFB code, the ground-state proper- 
ties of even-even ioo-i44gj^ isotopes such as binding energies per nucleon, two- 
neutron separation energies and nuclear charge radii have been calculated. The 
results are founded as to be successful. 
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